J /tp production in Indium-Indium collisions at 158 GeV/nucleon 
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The NA60 experiment studies muon pair production at the CERN SPS. In this letter we report on 
a precision measurement of 3/t[> in In-In collisions. We have studied the 3/t[> centrality distribution, 
and we have compared it with the one expected if absorption in cold nuclear matter were the only 
active suppression mechanism. For collisions involving more than ~80 participant nucleons, we 
find that an extra suppression is present. This result is in qualitative agreement with previous 
Pb-Pb measurements by the NA50 experiment, but no theoretical explanation is presently able to 
coherently describe both results. 
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The suppression of the charmonium states by color 
screening has been predicted as a signature of the phase 
transition from hadronic matter towards a Quark-Gluon 
Plasma [1]. Recently, it has been pointed out that 
the relatively loosely bound states ip' and \ c should in- 
deed melt for temperatures around the critical temper- 
ature T c , while the tightly bound 3/ip could survive, 
even if with strong medium- induced modifications, up to 
~ 2T C @, [1, [3] . However, the description of the evolution 
of a cc pair formed by gluon fusion at early times in the 
history of the collision, which may eventually lead to the 
formation of a bound state, is still not fully explained by 
theory. The influence of the medium, a percolatingpar- 
tonic condensate or a fully thermalized QGP p, or 
even a dense hadronic gas 0, H[ has been investigated in 
detail, but this physics topic is still largely data driven, 
and accurate experimental data are clearly needed. 

At the CERN SPS, the NA38 and NA50 experiments 
have already studied J/yj production in various collid- 
ing systems, including p-A |, S-U and Pb-Pb [u|. 
Proton-nucleus data provide an important reference, de- 
scribing the expected absorption of J /if) in cold nuclear 
matter. By comparing the centrality dependence of 
the J/ip yield observed in nucleus-nucleus collisions to 
this reference, one can look for suppression mechanisms 
connected with the formation of a strongly interacting 
medium. In particular, NA50 has observed, in Pb-Pb 



collisions, that below a certain centrality threshold the 
J /ip production is well described invoking nuclear absorp- 
tion as the only suppression mechanism; on the contrary, 
above that threshold, an extra suppression (also known 
as "anomalous" suppression) sets in. Such an interesting 
observation needs to be complemented by further sets of 
accurate measurements obtained with different collision 
systems. In this way one can study in more detail the on- 
set of the anomalous suppression, and understand which 
is the physics mechanism at its origin. 

The NA60 experiment has studied 3/yj production in 
In-In collisions at 158 GeV/nucleon at the CERN SPS, 
through its decay into two muons. The experimental ap- 
paratus is based on a muon spectrometer (MS), inherited 
from NA50 fl2j , used for triggering on muon pair produc- 
tion and for tracking purposes. A 12A.; thick hadron ab- 
sorber, mostly made of graphite, separates the MS from 
the target region, equipped with a beam tracker (BT) 
and a vertex tracker (VT), placed inside a 2.5 T dipole 
magnet. Finally, a Zero Degree Calorimeter (ZDC) 
provides an estimate of the centrality of the collisions. A 
more detailed description of the apparatus can be found 
in 14 , [l5| . The VT is a radiation-tolerant Si pixel de- 
tector, which tracks the charged particles produced in 
the collision (dV c h/df7 ~ 200 at midrapidity for central 
In-In interactions). By matching the tracks measured in 
the MS with the corresponding tracks in the VT, one 
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can access the kinematical variables of the muons before 
their distortion due to multiple scatterin g an d energy loss 
fluctuations in the hadron absorber 1(| 17j |. 

The results presented in this letter refer to the full In-In 
data sample collected by NA60, corresponding to ^230 
million events, taken at a beam intensity of 10 7 s _1 . Only 
the events having at least one interaction vertex with >4 
tracks attached to it are selected for the analysis. The 
distribution of the longitudinal coordinate of the vertices 
allows us to identify the target where the interaction took 
place with ~200 /im accuracy. We require the interaction 
vertex to lie, within this tolerance, in one of the 7 In 
sub-targets, 1.5 mm thick each and spaced by 7.5 mm. 
Matching of tracks in the VT with the MS tracks is then 
carried out in coordinate and momentum space. The 
matching efficiency, of the order of ^60% for events where 
a J/ip has been detected in the MS, shows no significant 
centrality dependence. In principle, MS tracks could be 
wrongly matched to VT tracks. However, in the J/ip 
mass region, this effect is negligible (<1%). 

Two kinds of event selection have been applied, corre- 
sponding to the two different analysis techniques detailed 
hereafter. In the first, the matching of the muon tracks is 
not performed, retaining in this way a larger event sam- 
ple. As a quality check, we perform a cut on the trans- 
verse distance between the extrapolated MS track at the 
target position and the beam axis, weighted by the mo- 
mentum of the track itself. The level of this cut has been 
set at 10% of the y 2 probability for each muon, a value 
that allows us to efficiently reject muons produced down- 
stream of the target. In the second selection, matching is 
applied and we require that the tracks matched to those 
detected in the MS originate from the reconstructed ver- 
tex or from the most upstream one, when more than one 
is found. In this way, we reject the small percentage 
(~4%) of events where the J/ip originates from a down- 
stream interaction of a nuclear fragment produced in a 
primary In-In collision. For both selections, we apply a 
beam pile-up rejection cut, based on the BT, requiring 
two subsequent ions to be separated in time by at least 
12 ns. In this way we avoid a superposition of signals 
in the read-out gate of the ZDC, that would bias the 
determination of the zero-degree energy -Ezdc • Finally, 
in order to reject events at the edges of the MS accep- 
tance we apply the kinematical cuts < y cm s < 1 and 
—0.5 < coshes < 0.5, where 6cs is the polar angle of 
the muons in the Collins-Soper reference frame. In the 
end, we are left with samples of ~45000 J/ip for the first 
event selection, and ^29000 for the second. 

The first analysis follows the approach used by the 
NA38/NA50 experiments. Namely, the J/ip yield is 
normalized to the corresponding Drell-Yan (DY), a 
hard process unaffected by final-state interactions in the 
medium [l8| . The ratio <jj / ody has the further ad- 
vantage of being free from systematic errors connected 
with the experimental inefficiencies and the integrated 



TABLE I: Values of o-j/^/ctdy, uncorrected for the J/ip decay 
branching ratio, as a function of centrality. 



Centrality bin 


(iVp art } 


"j/if/^DY 


Bzdc > 11 TeV 


63 


26.8±3.2 


7 < £zdc < 11 TeV 


123 


16.1±1.6 


Szdc < 7 TeV 


175 


17.8±1.6 



luminosity. It is extracted from a fit of the in- 
variant mass spectrum (in the region m MM > 2.2 GeV) 
to a superposition of the expected contributions, namely 
the J/ip and the ip' resonances, a continuum composed 
of Drell-Yan events and semi-muonic decays of D and D 
mesons, and a combinatorial background from tt and K 
decays. The expected mass shapes of the signals and 
their acceptances are evaluated through a Monte-Carlo 
simulation based on PYTHIA with GRV94LO [H 
parton distribution functions. The combinatorial back- 
ground has been estimated from the measured sample of 
like-sign pairs and its contribution is negligible. Fig. [1] 
shows the invariant mass spectrum, together with 

the result of the fit. 
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FIG. 1: The fit to the invariant mass spectrum, inte- 

grated over centrality. 

The ratios of cross sections, ctj/^/ctdyj uncorrected for 
the J/ip decay branching ratio, are given in Table U for 
three centrality bins corresponding to different values of 
(■^Vpart), the average number of participant nucleons in 
the bin. (N part ) has been obtained from £zdc using 
the Glauber model and taking into account the smearing 
induced by the detector resolution. The Drell-Yan yield 
refers to the mass interval 2.9 < m w < 4.5 GeV/c 2 . Its 
low statistics (320 events for > 4.2 GeV/c 2 ) limits 
the statistical significance of our result and prevents a 
finer binning in centrality. 

The values shown in Table Q] indicate that for semi- 
central and central collisions the crj /ctdy ratio is sig- 
nificantly lower than for peripheral reactions. In order to 
understand how much of the observed reduction is due 
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to cold nuclear matter effects, we have calculated, in the 
frame of the Glauber model, ctj m / ctdy as a function of 
centrality in a pure nuclear absorption scenario. This 
calculation requires as inputs the 3/tp absorption cross 
section in cold nuclear matter as well as the ratio of the 
J/ip and DY elementary production cross sections at 158 
GeV. Measurements performed by NA50 ll|, [2l| in p-A 
collisions at 450 GeV provide CTjfe=4.18 ± 0.35 mb and 
( cr j/i/j/ cr DY)45o=57.5 ± 0.8. The latter quantity has been 
rescaled to 158 GeV/nucleon incident energy, using a pro- 
cedure detailed in Ref. [ll|, obtaining (cj/v/ cr DY)i58 = 
35.7±3.0. Fig. [2] shows the measured <tj/^/<Jdy, divided 
by the calculated reference and plotted as a function of 
N paX f Even if statistical errors are large, a suppression 
signal of the J/ip beyond nuclear absorption can be seen. 




FIG. 2: Centrality dependence of the J/ip suppression mea- 
sured in In-In collisions. The stars correspond to the ratio 
between measured and expected ctj/^/ctdy, while the circles 
refer to the ratio between the measured J/ip yield and nu- 
clear absorption calculations. The box on the left shows the 
common systematic error, while the boxes around the points 
represent the error related to the centrality determination. 
See text for details. 



In the second analysis, we directly compare the mea- 
sured 3/tp yield to the centrality distribution of 3/tp cal- 
culated for the case of pure nuclear absorption, using the 
input parameters detailed above. In this approach, the 
measured J ftp yield has been obtained, in 1 TeV -Ezdc 
bins, by means of a simple fitting procedure that allows 
us to subtract the small amount of Drell-Yan (<4%) and 
combinatorial background (<1%) under the resonance 
peak. In Fig. [3] we compare the J/"0 distribution with 
the expectation from nuclear absorption. The relative 
normalization between data and the reference curve is 
not determined a priori; therefore we fix the ratio be- 
tween data and the nuclear absorption curve, integrated 
over centrality, to the same value (0.87 ± 0.05) as ob- 
tained within the previous analysis. It must be noted 
that events where a heavy nuclear fragment reinteracts in 
a downstream target have a smaller Ezuc value, since the 



fragment reinteraction removes some nucleons that would 
have otherwise reached the ZDC. The measured J/ip cen- 
trality distribution has been corrected for the small bias 
(<2%) induced by this effect. Finally, because of small 
inefficiencies of the BT, our data sample could be con- 
taminated by high-i?zDC events, where a non-interacting 
pile-up ion is superimposed to the interacting one. A 
Monte-Carlo simulation shows that this effect is negligi- 
ble in our analysis domain, i.e. Ezuc < 16 TeV (corre- 
sponding to Vp art > 50). 




FIG. 3: The J/ip -Ezdc distribution (circles), compared with 
expectations from nuclear absorption (line). 

Fig. [5] shows the ratio between the measured and ex- 
pected J ftp yield, rebinned in order to further reduce 
statistical fluctuations. Since statistical errors are very 
small (<2%), this analysis requires a careful estimate of 
systematic errors. We find that they are connected with 
the determination of the shape and normalization of the 
nuclear absorption reference, and with the calculation 
of Vpart starting from Ezbc- In particular, the uncer- 
tainties on (cj /^/ctdy)!?^ and af^^ give an 8% and 4% 
systematic error on the normalization of the absorption 
curve, respectively. We then have a 6% error, originating 
from the centrality integrated value of exj m / ctdy used for 
the normalization. Concerning centrality determination, 
by varying within errors the input parameters used in the 
Glauber model, we get a negligible influence on the nu- 
clear absorption reference, except for very central events 
(i?ZDC<3 TeV), where the size of the effect is —12%. 
Furthermore, the ZDC does not measure only spectator 
nucleons, but also a small amount of energy released by 
forward secondary particles emitted in the acceptance of 
the calorimeter (rj > 6.3). This contribution, important 
only for central collisions, is taken into account when cal- 
culating Vpart from Ezbc ■ By conservatively assuming a 
10% uncertainty on this quantity we get, for events with 
-Ezdc<3 TeV, a 9% error on the absorption curve. For 
more peripheral events the effect is negligible. Combin- 
ing the various sources in quadrature, we end up with 
a —11% systematic error, independent of centrality. On 
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top of that, the most central bins are affected by a fur- 
ther, sizeable systematic error relatively to the others. It 
must be noted that the systematic error plotted in Fig. [21 
except for the fraction due to the 6% normalization error 
quoted above, also affects the results on oj/^/ctdy- 

The results obtained by NA60 in In-In collisions show 
that in the region 50 < N pait < 100 an anomalous sup- 
pression of the J/ip sets in. Taking into account the iV par t 
smearing due to the -Ezdc resolution, the effect seen is 
compatible with the occurrence of a ~15% drop of the 
J/ip yield at iVp ar t ~ 80, followed by a more or less con- 
stant behaviour. When expressed in terms of the Bjorken 
energy density, the onset of the anomalous suppression 
roughly corresponds to 1.5 GeV/fm 3 (using t =1 fm/c 
and the VENUS [22| estimate for the charged multiplic- 
ity as a function of centrality). In Fig. 3] we compare 
our result with the J/ip suppression pattern obtained by 
NA50 in Pb-Pb collisions n|. The systematic errors 
on the determination of the nuclear absorption reference 
from the p-A data sample amount to ±9% and are not 
shown in this comparison plot since they affect Pb-Pb 
and In-In results in a similar way. Within errors, the two 
patterns look compatible in the iV part region explored by 
both systems, indicating that -/V part might be, at SPS en- 
ergy, a scaling variable for the anomalous suppression. 
A detailed investigation of the scaling properties of J/ip 
suppression as a function of several centrality variables 
would give valuable insights into the origin of the ob- 
served effect. However, a meaningful comparison would 
require Pb-Pb results with error bars similar to the ones 
obtained for the In-In analysis. 
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FIG. 4: Comparison between the In-In (NA60, circles) and 
Pb-Pb (NA50, triangles) suppression patterns. The box on 
the left shows the 6% systematic error related to the normal- 
ization procedure of the In-In points. 

Several theoretical predictions, tuned on already avail- 
able Pb-Pb results from NA50, were formulated for J/ip 
suppression in In-In collisions 0, 0, Hil- We find that 
none of them is able to quantitatively rep roduce the sup- 
pression pattern measured by NA60_|24j. Recent results 
from the PHENIX Collaboration [25J have shown that 



also in Au-Au collisions at -^5=200 GeV/nucleon the J/ip 
is suppressed beyond nuclear absorption, and that the 
suppression is larger at forward rapidities. A coherent 
interpretation of the results at SPS and RHIC energies is 
now mandatory in order to understand the physics mech- 
anisms affecting charmonia in a dense partonic/hadronic 
environment. The results obtained by NA60 represent 
the most accurate measurement of J/ip suppression in 
nuclear collisions available today and are a key element 
to strictly constrain theoretical models. Further studies 
on the y and px dependence of the J/ip suppression are 
underway and will be the subject of future publications. 

In summary, we have measured J/ip suppression in 
In-In collisions at 158 GeV/nucleon. Comparing the J/ip 
centrality distribution with the expectation from a pure 
nuclear absorption scenario, we find an anomalous sup- 
pression that sets in for A part ~ 80, and saturates for 
more central events. The statistical errors are negligi- 
ble (of the order of 2%). Most of the systematic errors 
are centrality independent and therefore do not affect the 
measured shape of the J/ip suppression pattern. None of 
the existing theoretical models, tuned on the measured 
J/tp suppression in Pb-Pb collisions, is able to quantita- 
tively reproduce the results shown in this letter. 
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Foundation, Lisbon and the Swiss Fund Kidagan. 
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